Entropic transport - A test bed for the 
Fick-Jacobs approximation 

By p. Sekhar Burada, Gerhard Schmid and Peter Hanggi 

Institut fiir Physik, Universitdt Augsburg, Universitdtsstr. 1, D-86135 Augsburg, 

Germany 

Biased diffusive transport of Brownian particles through irregularly shaped, narrow 
confining quasi-one-dimensional structures is investigated. The complexity of the 
higher dimensional diffusive dynamics is reduced by means of the so-called Fick- 
Jacobs approximation, yielding an effective one-dimensional stochastic dynamics. 
Accordingly, the elimination of transverse, equilibrated degrees of freedom stem- 
ming from geometrical confinements and/or bottlenecks cause entropic potential 
barriers which the particles have to overcome when moving forward noisily. The ap- 
plicability and the validity of the reduced kinetic description is tested by comparing 
the approximation with Brownian dynamics simulations in full configuration space. 
This non-equilibrium transport in such quasi-one-dimensional irregular structures 
implies for moderate-to-strong bias a characteristic violation of the Sutherland- 
Einstein fluctuation-dissipation relation. 
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1. Introduction 

Diffusion of Brownian particles through narrow, tortuous confining structures such 
as micro- and nano-pores, zeolites, biological cells and microfluidic devices plays a 
prominent role in the dynamical characterization of these systems (Barrer, 1978; 
Berezhkovskii & Bezrukov, 2005; Hille, 2001; Kettner et al, 2000; Liu et al, 1999; 
Matthias & Miiller, 2003; Miiller et al, 2000; Nixon & Slater, 2002; Siwy et al, 
2005; Volkmuth & Austin, 1992). Effective control schemes for transport in these 
systems requires a detailed understanding of the diffusive mechanisms involving 
small objects and, in this regard, an operative measure to gauge the role of fluc- 
tuations. The study of these transport phenomena is in many respects equivalent 
to an investigation of geometrically constrained Brownian dynamics (Burada et al, 
2009; Hanggi & Marchesoni, 2009; Mazo, 2002). With this work we focus on the 
stochastic transport of small sized particles in confined geometries and the feasi- 
bility of the so-called Fick-Jacobs (FJ) approximation to describe the steady-state 
particle densities. Restricting the volume of the configuration space available to the 
diffusing particles by means of confining boundaries or obstacles discloses intriguing 
entropic phenomena (Liu et al, 1999). 

The driven transport of charged particles across bottlenecks (Burada et al, 
2009), such as ion transport through artificial nanopores or artificial ion pumps 
(Kosinska et al, 2008; Siwy et al, 2005) or in biological channels (Berezhkovskii & 
Bezrukov, 2005) are more familiar systems where diffusive transport is regulated by 
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entropic barriers. Similarly, the operation of artificial Brownian motors and molec- 
ular machines relies as well on a mutual interplay amongst diffusion and binding 
action by energetic or, more relevant in the present context, entropic barriers (As- 
tumian & Hanggi, 2002; Burada et al., 2009; Derenyi & Astumian, 1998; Hanggi & 
Marchesoni, 2009; Reimann & Hanggi, 2002). 

The outline of this work is as follows: In Sec. 2 we introduce our model and 
formulate the mathematical formalism needed to model the diffusion of a Brownian 
particle immersed in a confined medium. In Sec. 3 we present the FJ approximation 
and compute the entropic effects on the particle transport and on the steady-state 
probability density in the presence of an applied force in transport direction. In 
Sec. 4 we compare the numerical precise 2D simulation results with those obtained 
from applying the FJ approximation. In Sec. 5 we discuss the effective lateral dif- 
fusion and test the Sutherland-Einstein fluctuation-dissipation relation. Sec. 6 pro- 
vides a discussion of our main findings. 

2. Overdamped system dynamics 

Generic mass transport through confined structures such as irregular pores and 
channels, c.f. the one depicted with Fig. 1, is governed by the transport of suspended 
Brownian particles subjected to an externally applied potential V{f}. Generally, 
the dynamics of the Brownian particle inside the medium can be well described 
by a Langevin dynamics in the over-damped limit (Purcell, 1977), with reflecting 
boundary conditions at the channel walls. The stochastic dynamics then reads 



(2.1) 



where r is the position vector of a Brownian particle at time t, q denotes the 
friction coefficient, fee is the Boltzmann constant and T refers to the environmental 
temperature. Thermal fluctuations due to the coupling of the Brownian particle 
to the environment are modeled by Gaussian white noise with zero mean and an 
auto-correlation function obeying the Sutherland-Einstein fluctuation-dissipation 
relation (Hanggi & Marchesoni, 2005): 



= 2<5,, 6{t-t') for i, j = x,y,z. 



(2.2) 




Figure 1. Schematic diagram of a channel confining the motion of laterally forced (with 
strength F)Brownian particles. The half-width a; is a periodic function of a; with periodicity 
L. 



For simplicity we consider the dynamics of a Brownian particle that is subjected 
to an constant force F = Fex acting along the direction of the channel axis (here 
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in a;-direction) . The Langevin equation for the over-damped dynamics then reads: 



T] =F+ V??fcBT|(f) , (2.3) 

with reflecting (i.e. no across flow) boundary conditions implied at the channels 
walls which confine the Brownian particles within the channel geometry. 

In order to further simplify the treatment of this set up we introduce dimension- 
less variables. We measure all lengths in units of the period length L, i.e., r = Lx, 
where x denotes the dimcnsionlcss position vector of the particle. As the unit of 
time r wc choose twice the time the particle takes to diffusively overcome the dis- 
tance L, which is given by t = L'^r]/{kBT), i.e. i = rt (Burada et al, 20086). In 
these dimensionless variables the Langevin dynamics assumes the form 

f = / + !'(*)• (2.4) 

where {^{t)) = 0, {^i{t) ^j{t')) = 2 Sij 5{t — t') for i,j = x,y, z, and the dimensionless 
force becomes 

f = fe, and/=^. (2.5) 

The performed dimensionless scaling parameter / characterizes the force as the ratio 
between the work LF done on the particle along a distance of the period length 
L and the thermal energy k^T . We anticipate here the fact that in the case of 
diffusion occurring in purely energetic potential landscapes the driving force F and 
the temperature T are independent variables; in contrast, in systems with entropic 
features these two quantities become coupled (Reguera et al. 2006). In order to 
adjust a certain value of / one can modify either the force strength F or adjust the 
noise intensity /cbT. 

The corresponding Fokker-Plank equation describing the time evolution of the 
probability density P{x,t) takes the form (Hanggi & Thomas, 1982; Risken, 1989): 

^^1^ = -V-J(f,i), (2.6a) 

where J{x,t) is the probability current: 

J(f,i)= (/-V)P(x,i). (2.6b) 

Note that, for channels with similar geometry which are related by a scale trans- 
formation r Xf, A > 0, the transport properties are determined by a single 
dimensionless parameter / which subsumes the respective period length, the exter- 
nal force and the temperature of the surrounding fluid. 

The no-flow condition beyond the channel walls leads to a vanishing probability 
current at those boundaries. Consequently, due to the impenetrability of the channel 
walls, the normal component of the probability current J{x,t) vanishes at those 
boundaries. Thus, the boundary conditions at the channel walls are given by 



J{x, t) ■ n = , X € channel wall , (2.7) 
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where n denotes the normal vector at the channel walls. 

The boundary of a 2D periodic channel, which is mirror symmetric about the 
a;-axis, is given by the dimensionless periodic functions y = ±u){x), i.e., lo{x + 1) = 
clj{x) for all x, where x and y are the cartesian components of x. In this case, the 
boundary condition reads 



da;(a;) 
dx 



fP{x,y,t)- 



dP{x,y,t) 
dx 



I dP{x,y,t) 
dy 



(2.8) 



at y = ±w(a;). Except for a straight channel with uj = const, there are no periodic 
channel shapes for which an exact analytical solution of the Fokker-Planck equation 
(2.6a, 2.6b) with the elaborate boimdary conditions in (2.8) is presently known. 
Approximate solutions though can be obtained on the basis of a one-dimensional 
diffusion problem proceeding in an effective potential. Narrow channel openings, 
which act as geometric hindrances in the original system, then manifest themselves 
as entropic barriers within an effective one-dimensional diffusive FJ approximation 
(Burada et al., 2007; Jacobs, 1967; Kalinay & Percus, 2006; Reguera et al, 2006; 
Reguera & Rubi, 2001; Zwanzig, 1992). 



3. Equilibration in transverse channel directions: the 
Fick-Jacobs approximation 

In the absence of an external force, i.e. for / = 0, it was shown (Jacobs, 1967; 
Kalinay & Percus, 2006; Reguera & Rubi, 2001; Zwanzig, 1992) that the dynamics 
of Brownian particles in confined structures (such as the one depicted in Fig.l) can 
be described approximatively by the FJ equation; i.e.. 



^P(x, t) = ^ D{x) e-^(-) ^e^(-)p(x, t) . (3.1) 
ot ox ox 

This ID equation is obtained from the full 2D Smoluchowski equation upon the 
elimination of the transverse y spatial coordinate degree of freedom by assum- 
ing a much faster equilibration in that channel direction than in the longitudinal 
one. An analogous reduction mechanism has been used for the transport of neu- 
trons through nuclear reactors (Beckurts & Wirtz, 1964). In the equation (3.1), 

P{x, t) = /"^^;|,) P{x, y, t) dy denotes the marginal probability density along the 
axis of the channel. A{x) corresponds to the potential of mean force which equals 
for the considered situation the free energy, i.e. A{x) = E{x)—S{x) = 0— lna;(a;). We 
note that for three dimensional channels an analogue approximate Fokker-Planck 
equation holds in which the function u^x) is to be replaced by 7rcL'^(.x) (i.e. the 
area of the corresponding cross-section). In the original work by Jacobs (Jacobs, 
1967) the ID diffusion coefficient D{x) is constant and equals the bare diffusion 
constant which assiimcs unity in present dimensionless variables. However, intro- 
ducing the x-dependent diffusion coefficient considerably improves the accuracy of 
the kinetic equation, extending its validity to more winding structures (Burada et 
al., 2007; Reguera & Rubi 2001; Zwanzig, 1992). The expression for D{x) reads (in 
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Figure 2. Sketch of the 2D channel and the effective one-dimensional potential: the Fick- 
-Jacobs (FJ) approximation allows for a reduction of the 2D Brownian dynamics within 
the periodic channel (periodicity: L) to an approximate ID Brownian dynamics with an 
effective potential which is given by the free energy function A{x). In the presence of an 
applied bias A{x) has the form of a tiled periodic potential with a barrier height of AA 
which depends on the temperature T. 



dimensionless units) 

D{x) = ^ ]. , (3.2) 

where a = 1/3, 1/2 for two and three dimensions, respectively, has been shown to 
appropriately account for the curvature effects of the confining walls (Burada et al, 
20086; Reguera & Rubi, 2001). uj'{x) indicates the first derivative of the boundary 
function uj{x) with respect to x. 

In the presence of a constant force F along the direction of the channel the FJ 
diffusion equation (3.1) can be recast into the form (Burada et ai, 2007; Burada et 
al, 20086, Reguera et at, 2006): 



with the dimensionless free energy A{x) := E{x) — S{x) = —fx — \iiuj(x). In terms 
of the original unsealed physical variables the energy is E = k-QTE{x) = —Fx {x — 
xL) and the dimensional entropic contribution reads 5* = kBTS{x) = k^T hiuj{x). 
For a periodic channel arrangement this free energy assumes the form of a tilted 
periodic potential, see Fig. 2. In the absence of a force the free energy is purely 
entropic and Eq. (3.3) reduces to the FJ equation (3.1). On the other hand, for 
a straight channel the entropic contribution vanishes and the particles are solely 
driven by the externally applied force. 

Remarkably, the temperature T dictates the strength of the effective potential. 
An increase in temperature causes an increase in barrier height while for purely 
energetic systems the barrier height is independent of the temperature (Hanggi et 
al, 1990). 
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(a) Steady- state probability density 

Formally, the steady-state density of the particles is obtained in the limit f — > oo, 
i.e. P^*(x) = limt^oo As a consequence, ^P^*(a;) = 0. An expression for 
the steady-state density can be derived from Eq. (3.1), using arguments detailed in 
the Appendix. Using the main result in Eq. (A 12) one obtains 

P^\x)= /^^'-^^ , (3.4) 

/ I{z,f)dz 
Jo 

where 

depends on the dimensionless position x, the force / and via the position dependent 
diffusion coefficient on the shape of the tube given in terms of the shape function 
u!{x) and its first derivative, cf. Eq. (3.2). Note, that the probability density P^^{x) 
is normalized on the unit interval. 



(6) Nonlinear mobility 

The primary quantity of particle transport through periodic channels is the 

average particle current, (x) , or cquivalently the nonlinear mobility, which is defined 
as the ratio between the average particle current and the applied force /. For the 
average particle current we derive an expression which is similar to the Stratonovich 
formula for the current occurring in titled periodic energy landscapes, but hero 
with a spatially dependent diffusion coefficient (Burada et al., 2007). A detailed 
derivation of this expression is given in the Appendix (Appendix A), cf. Eq. (All). 
Hence, we obtain the nonlinear mobility for a 2D or 3D channel: 



M/) = ^ = ) . (3.6) 

^ ^ / I{zJ)dz 
Jo 

with I{z, /) given in Eq. (3.5). 

4. Precise numerics for a 2D channel geometry 

The steady-state density and the average particle current, predicted analytically 
above, has been compared with Brownian dynamic simulations performed by a 
numerical integration of the Langevin equation Eq. (2.4), using the stochastic Euler- 
algorithm. The shape of the exemplarily taken 2D channel is described by 

a;(a;) := asin(27ra;) -I- 6, (4.1) 

where b > a. The sum and difference of the two parameters a + b and b — a yield 
half of the maximal and the minimal width of the channel, respectively. Moreover, a 
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Figure 3. The normalized steady-state probability density of particles along the propaga- 
tion direction is depicted for different force values / (see in panel (b)) for two different 
2D-channels (see inset in panel (a)), with the scaled half- width shape function given by 
uo{x) — a [sin(27r(a; — 0.25)) + 1.02] (the shift —0.25 ensures that the bottlenecks are lo- 
cated at and 1). For a — 1 the maximal and minimal channel widths are 4.04 and 0.04, 
respectively; likewise for a = l/27r they are 6.43 ■ 10~^ and 6.37 ■ 10"'^, respectively. The 
Solid lines correspond to the steady-state probability density obtained from ID Fick- Jacob 
approximation, Eqs. (3.4, 3.5), and symbols correspond to 2D numerical simulations, see 
Eq. (4.2). 



controls the slope of the channel walls which in turn determines the one-dimensional 
diffusion coefficient D{x). 

For the considered channel configuration, cf. Eq. (4.1), the boundary condition 
becomes uj{x) = a [sin(27rx) + k], where k = b/a = 1.02 throughout this paper. 
For a we chose values between 1 and l/27r. In all cases the width of the widest 
opening within the channel is larger by a factor of 100 than the width at the 
narrowest opening. One may therefore do expect rather strong entropic effects for 
these channel geometries. 



(a) Stationary probability densities 

We have evaluated the stationary probability density P^^{x,y), in the long time 
limit, by mapping all particle positions onto the primitive cell by translation in 
longitudinal channel direction. Consequently, da; /"f,^^-, dy P^*^{x,y) — 1. Fig. 3 
(solid lines) depicts the normalized steady-state probability density in x-direction, 
for various scaling parameter values, derived from the reduced one-dimensional FJ 
result, Eq. (A 12), and is compared with the numerical simulations for the exact 
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(a) 

a = 1.0 




f = 0.2 
/ = 3.0 
/ = 10.0 
/ = 50.0 



(b) 

a = \/2tt 






0.75 
0.5 
0.25 




y 

Figure 4. The normalized steady-state probability density of particles in t/-direction P^^{y), 
cf. Eq. (4.3), taken at the a::-position of maximal channel width, is depicted for different 
values of the scaling parameter / (see in panel (b)) and different channel geometries (see 
inset in panel (a)) for the boundary function uij{x) = a [sin(27r(a; — 0.25)) + 1.02] (the 
shift —0.25 ensures that the bottlenecks are located at and 1). The maximal channel 
widths for the two structures are: 2!i;inax ~ 4.04 for a = 1 and 2a;niax = 6.43 ■ 10~^ for 
a = 1/2-n. The symbols correspond to 2D numerical simulations, see Eq. (4.3). For large 
scaled force values / the numerically obtained steady-state densities deviate from the 
uniform, i.e. box-like, density (shown by the black solid line). This indicates the failure of 
the equilibration assumption on which the Fick-Jacobs-approximation relies. 



expression: 



P''*(x) 



P'\x,y) dy 



dx 



(4.2) 



P''ix,y) dy 



Note, that steady-state marginal density P^'(x) is normalized on the primitive ceU. 

At small scaling parameter values / the ID steady-state density given by Eq. (3.4) 
is in very good agreement with those obtained from numerical simulations, see 
Fig. 3. This holds true for rather arbitrary channel geometry (not shown). How- 
ever, the comparison fails for large /-values of the scaling parameter or for more 
winding structures corresponding to larger a- values. When increasing the force, the 
maximum of P^*^{x) is shifted towards the exit of the cell and the particles mostly 
accumulate in front of the bottleneck, see Fig. 3(a), and the ID kinetic description 
starts to fail in that forward bottleneck x-region. However, by decreasing a of the ge- 
ometric channel shape function the accuracy of the FJ approximation considerably 
improves up to very large force values /, see Fig. 3(b). 
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As a common feature one observe that for the two chosen geometric structures 
that in the large force regime the numerically obtained P^^{x) is essentially con- 
stant over a wide range of x-valucs, indicating a minor influence of the shape of 
the structure on the dynamics of the laterally forward-forced particles. In this sit- 
uation, the thermal noise plays a minor role and the deterministic dynamics (with 
diffusion set to zero) of the diffusive equation provides a good starting point. Put 
differently, at strong longitudinal driving strength the correction in the diffusion 
coefficient leading to a spatial dependency, i.e. D{x), overestimates the role of the 
entropic effects and consequently the FJ approximation starts failing over extended 
a;-regimes. 

The reasons for the failure of the FJ approximation for large forces becomes ob- 
vious when checking the equilibration assumption in transverse channel direction. 
From our simulations, we can actually analyze the validity of the hypotheses of equi- 
libration in the transverse direction on which the FJ description relies. A detailed 
analysis is provided by testing the normalized steady-state probability density in 
the transverse direction at a given x-position, i.e. 

P^'iy)--= jy''^ dy. (4.3) 



/u>(x) 



In Fig. 4, we depict the steady-state probability density at the position of max- 
imal channel width. For small values of the scaling parameter /, the P^{y) is very 
flat, indicating an almost ideal homogeneous equilibration in the transverse direc- 
tion, as required by the FJ approximation scheme. However, at large force strengths 
/ the Brownian particles concentrate along the axis of the channel with y = Q. In 
this situation, the assumption of equilibration along the transverse direction fails 
and the density peaks around the y = value. The particles can only feel the pres- 
ence of the boundaries when they are close to the bottlenecks. Hence, in the limit of 
very large force values, the influence of the entropic barriers practically disappears. 

(6) Nonlinear mobility 

The average particle current was derived from an ensemble-average using 3 ■ 10^ 
trajectories: 

{i) = lim ^ . (4.4) 

Fig. 5 shows the nonlinear mobility as a function of the scaling parameter /. We 
note that the transport in one-dimensional periodic energetic potentials distinctly 
differs from the one occurring in one-dimensional periodic systems in presence of 
entropic barriers (Reguera et ai, 2006). The fimdamental difference lies in the 
temperature dependence of these barrier shapes. Decreasing the temperature in 
an energetic periodic potential decreases the transition rates from one cell to the 
neighboring one by decreasing the Arrhenius factor exp{— Ay/(fcBT')} where AV^ 
denotes the activation energy necessary to proceed a period (Hanggi et al., 1990). 
Hence decreasing the temperature yields a decreasing nonlinear mobility. For a 
one-dimensional periodic system with an entropic free energy (or entropic potential 
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1 10 100 

/ := FL/kBT 

Figure 5. The numerically simulated (symbols) and analytically calculated (cf. Eq. (3.6) 
- lines) dependence of the scaled nonlinear mobility n{f) vs. the dimensionless force 
/ = FL/ksT is depicted for two 2D channel geometries. For both channels the scaled 
half-width is given by u}{x) — a [sin(27ra;) + 1.02]; a = 1 : circles and solid line, 
a = 1/{2-k) : diamonds and dashed line. The dotted line indicates the deterministic limit 
M(/) = {i)/f = 1- 

of mean force), a decrease of temperature results, however, in an increase of the 
dimensionless force parameter / and consequently in a monotonic increase of the 
nonlinear mobility, cf. Fig. 5. 

The dependence of the dynamics on the geometry parameter a nicely reflects the 
entropic effects on the mobility: A channel with a larger a- value has wider openings 
and thus provides more configuration space where the particle can sojourn. This 
longer residence time within a period of the channel diminishes the throughput and 
consequently the mobility. This is corroborated by the results of our calculations 
depicted in Fig. 5. For all values of /, an increase in value of a leads to a decrease 
in the mobility. This holds not only in regimes for which the FJ equation is valid, 
but also for large values of / where the approximation fails. For very large values 
of the scaling parameter / the nonlinear mobility approaches the value 1, i.e., it 
agrees with the deterministic strong driving limit. 

By means of the nonlinear mobility a detailed comparison between 2D simulation 
results and the analytic results, cf. Eq. (3.6), enables one to determine validity 
criteria for the FJ approximation, for further details see in Ref. (Burada et ai, 
2007; Burada et al, 20086). 

5. Effective diffusion and the Sutherland-Einstein relation 

A validity of a nonlinear Sutherland-Einstein relation implies that in physical units 
we can relate the nonlinear mobility IJ,{F) directly to the nonlinear, effective x- 
diffusion Deff{F), reading 

D.ifiF) = ^IiF) kuT , (5.1) 

Put differently, the effective diffusion coefficient D^g for the diflFusive spreading 
along the longitudinal channel direction would then solely be determined by the 
nonlinear mobility discussed above and the environmental temperature T. 
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Figure 6. The ratio of the effective Diffusion Dgg and nonlinear mobility times the 
thermal energy UbT is depicted as a function of the scaling parameter / for two channel 
geometries: ij{x) = a[sin(27rx) + 1.02] with a = 1 (circles) and a = l/27r (diamonds). The 
dashed horizontal line at 1 indicates the validity of the Sutherland-Einstein relation in 
this non-equilibrium situation: DcB /{fiksT) = 1. A deviation from this line consequently 
marks the breakdown of this relationship. The inset depicts the ratio of the effective 
diffusion Des and bulk diffusion constant, being kaT/ri. A ratio larger than 1 (dashed 
horizontal line) indicates a characteristic enhancement of the effective a;-diffusion. 

A validity of this relation would then imply a monotonic increase towards the 
entropic-free diffusion limit, i.e. D^ff = ksT/r]. The latter is being approached in 
the strong forcing hmit where entropic effects cease to play a significant role. Such 
a monotonic behavior, however, is not observed from the numerical simulations for 

the effective .T-diffusion coefficient (Burada et ai, 20086; Regucra et ai, 2006). It is 
defined as the ratio between the asymptotic behavior of the variance of the position 
variable and the elapsed time t; i.e., 

D^s = hm ^ ^ ^' ■ (5.2) 

Interestingly, the dependence of the effective diffusion coefficient on the scaling 
parameter exhibits a bell-shaped behavior, cf. inset of Fig. 6, thus indicating a fail- 
ure of the Sutherland-Einstein relation in this moderate-to-strong driving regime. 
This break down of the Sutherland-Einstein relation can also be detected within 
the FJ description (not shown in Fig. 6): The FJ approximation for this effective 
x-diffusion as well yields a non-monotonic dependence of the effective a;-diffusion 
coefficient on the scaling parameter / exhibiting a peak value exceeding the bulk 
diffusion coefficient Dq = rj/ksT (Burada et ai, 20086; Reguera et ai, 2006). 

For a detailed comparison, we depict the ratio of numerically obtained DeS 
and {iik^T) in Fig. 6. Surprisingly, it turns out that such a Sutherland-Einstein 
relation, Eq. (5.1), holds true in terms of the effective mobility in the small forcing 
limit / — > ; i.e. in the linear response regime. It increasingly fails, however, for 
increasing bias strength F. At very strong bias, i.e., / — * oo the biased diffusion 
becomes effectively " free" from entropic effects and cxpectedly approaches the free 
limit, given by ksT/r], which renders the original, linear Sutherland-Einstein result 
in terms of the F-independent mobility /i = l/jj. Put differently, the influence of 
entropic barriers caused by the bottlenecks becomes negligible at strong bias. Vice 
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versa, the bell-shaped behavior of the ratio depicted with Fig. 6 reflects the fact 
that this effective diffusion is not increasing monotonically but rather exhibits an 
enhancement of effective diffusion at moderate bias (or scaling) values /, c.f. in the 
inset of Fig. 6. 



6. Conclusions 

In summary, we demonstrated the applicability of the equilibration approxima- 
tion in describing biased diffusive transport occurring in narrow, irregularly shaped 
one-dimensional channel structures. The Fick-Jacobs (FJ) description which relies 
on the equilibration assumption allows for a treatment of the; dynamics within an 
effective one-dimensional kinetic equation of the Smoluchowski form. Bottlenecks 
and other confining restrictions of available configuration space yield within this 
approximation an effective ID diffusion equation exhibiting entropic barriers. Due 
to the intrinsic temperature dependence of the underlying entropic free energy con- 
tribution one finds for the transport phenomena in periodic channels possessing 
varying cross sections features that arc radically different from conventional trans- 
port occurring in energetic periodic potential landscapes. 

The most striking difference between these two physical situations is that for a 
fixed channel geometry the dynamics is characterized by a single scaling parame- 
ter / = FL/{k^T) which combines the external force F causing a drift, the period 
length L of the channel, and the thermal energy k-oT. The latter presents a measure 
of the strength of the acting fluctuating thermal forces. This leads to an opposite 
temperature dependence of the mobility: While the mobility of a particle in an 
energetic potential increases with increasing temperature the mobility of a particle 
undergoing biased diffusion in an irregular channel decreases. The incorporation of 
the spatial variation of the channel width in terms of an entropic free energy contri- 
bution allows for a quantitative understanding of the dependence of the transport 
properties, like the nonlinear mobility, on parameters like force strength, channel 
topology or temperature. Moreover, the lateral steady-state probability densities 
P^^{x) can be evaluated in analytical closed form within the reduced kinetic FJ 
approximation, see in the Appendix. 

Such an effective one-dimensional reduction of a complex diffusion dynamics 
with intricate boundary conditions at the confining walls certainly proves useful and 
beneficial for the quantitative description, design and control of diffusive transport 
along tortuous pores and alikc\ The latter situation dictates the stochastic far- 
from-equilibrium transport in a great variety of biological and structured synthetic 
pores and confining cavities, such as buckyballs, zeolites and alike. As an example, 
this FJ approximation has successfully been used in describing the phenomenon 
of Stochastic Resonance (Hanggi, 2002; Gammaitoni et al., 1998) in a 2D system 
exhibiting an entropic barrier (Burada et al, 2008a). 
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Appendix A. Steady-state Current and Probability Density 

In this appendix we derive the steady-state solution for the effective, one-dimensional 
dimensionless FJ (Smoluchowski-type) equation, Eq. (3.1), 

^P(x, t)^^ D{x) e-^(-) Ae^(-)p(x, t) , (A 1) 

where A{x) denotes the free energy function, A{x) = — f x — \nu!{x), withti;(x+l) = 
Lo{x). Eq. (A 1) results from the probability continuity equation; 

^P{x,t) = -^J{x,t), (A 2) 

with the probability current J{x, t) reading: 

J{x, t)= - D{x) e-^(") ^e^(^^P(a;, t) . (A 3) 

In the case of a tilted periodic potential, i.e. A{x + 1) = A{x) — /, it is convenient 
to define the reduced probability density and the corresponding reduced current; 
i.e., 

P(x,i)=^P(n + a;,t), (A 4) 

n 

J{x, t) = ^J{n + x,t), nGZ. (A 5) 

n 

By definition these functions are periodic with periodicity i = 1, P(.x + 1, i) = 
P{x, t), J{x+1, t) = J{x,t). The same holds true for the spatial dependent diffusion 
coefficient, i.e. D{x + 1) = D{x). Moreover, P{x,t) and J{x,t) obey the continuity 
equation, Eq. (A 2), and P{x,t) is normalized on an interval {x,x + 1), provided 
that P(x, t) is normalized, i.e. P{x, t) dx = 1. 

In the steady-state limit the probability current assumes a constant, i.e., J{x, t) —>■ 
J. Thus Eq. (A 3) becomes 

J= -P)(x)e-^(^)^e^(^)P'''(a;)- (A6) 
ox 

Multiplying both sides of Eq. (A 6) by e"^^''^ /D{x), and integrating over a period 
i = 1 we obtain 

which simplifies to 

rx+l A{x') 

J ^^dx' = P^\x){l - e-Oe^(-). (A8) 

Upon rearranging the terms on the right hand side and integrating once more over 
a period, i.e. from to 1, we find the first result 

J=—^ J /^ • (A9) 



h Jx D{x') 
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Hereby, we made use of the normalization condition of the stationary probability, 
i.e., P^^{x)dx = 1. The general relation between the steady-state probability 
current and the steady-state average particle current {{x)) is 

(i) = / J da; (A 10) 

Jo 

which implies that {x) = J. Thus, the transport current is given by the first main 
result, reading: 

(a;) = J . (A 11) 

^ ' /■! ^K+l A(x) ^ ' 



By substituting Eq. (A 9) back into the Eq. (A 8) we obtain for the steady-state 
probability density in a;-direction the second main result: 



/•rr + l A{x') 

/ D(x') 

P'\^) = ■ (A 12) 



r 

J X 

fl fX+1 A{x') 

Jo Jx D{x') 
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